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Abstract 



J2L|| We investigate complex-temperature singularities in the Ising model on the triangular lattice. 

pH I Extending an earlier analysis of the low-temperature series expansions for the (zero-field) 

susceptibility x by Guttmann |]I| to include the use of differential approximants, we obtain 
K> , further evidence in support of his conclusion that the exponent describing the divergence in 

^ . X at u = Me = — 1/3 (where u = e~^^) is 7g = 5/4 and refine his estimate of the critical 

amplitude. We discuss the remarkable nature of this singularity, at which the spontaneous 
magnetisation diverges (with exponent (3e = —1/8) and show that it lies at the endpoint 
of a singular line segment constituting part of the natural boundaries of the free energy in 
the complex u plane. Using exact results, we find that the specific heat has a divergent 
singularity at u = —1/3 with exponent a'^ = 1, so that the relation a'^ + 2j3e + 7g = 2 is 
satisfied. We also study the singularity at u = Ug = —1, where M vanishes (with Ps = 3/8) 
and C diverges logarithmically (with a^ = a^ = 0). 
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1 Introduction 

Although exact closed-form expressions for the (zero-field) free energy and spontaneous mag- 
netisation of the two-dimensional Ising model have been calculated, no such expression is 
known for the susceptibility. One can have at least two views on this. The first is that since 
the problem has been unsolved for 50 years, it is too difficult and one should consider other 
problems more amenable to solution. A second view is that some progress on even longstand- 
ing problems can be made with the input of new information. In support of this view, one 
could cite, for example, the recent progress on proving Fermat's last theorem, which has 
been unsolved for more than 350 years since its first assertion in 1637. We shall adopt the 
second view here. In this context, every additional piece of information on the susceptibility 
is of value, especially insofar as it specifies properties which an exact solution must satisfy. In 
particular, it is of interest to study the susceptibility as an analytic function in the complex 
temperature plane. Complex-temperature singularities of the susceptibility have been known 
for many years. When early low-temperature series expansions were calculated for the Ising 
model by the King's College group [0], it was found that for the simple cubic, body-centered 
cubic, and face-centered cubic lattices, these have complex-temperature singularities closer 
to the origin than the position of the physical critical point, which precluded the use of a 
ratio test to determine the location of this physical critical point 0, ^ . For the triangular 
lattice, it was found that the series expansions showed a complex-temperature singularity at 

the point [|, || 

1 , , 

U = Ue = -- (1.1) 

on the opposite side of the origin, and an equal distance from it, compared to the physical 
critical point Mc = 1/3 for this lattice, where 

u = e-^^ (1.2) 

and K denotes the spin-spin coupling (see below for notation). Transformations of the series 
variables to map the unphysical singularities farther away from the origin than the physical 
critical point and methods of analysis such as dlog Fade approximants enabled one to obtain 
useful information on critical behaviour from low-temperature series (for an early review, see 
Ref. 0). In particular, as part of a general analysis of the susceptibility of the Ising model 
on 2D lattices |jl| using high- and low-temperature series, Guttmann studied the complex- 
temperature singularity at m = — 1/3 for the triangular lattice. To obtain a series with good 
sensitivity to this singularity, he subtracted the dominant physical singularity and its leading 
confiuent correction, and analysed the resultant series using ratio and dlog Fade methods. 



From this analysis, Guttmann concluded that as u approaches Ue from the direction of the 
origin, the susceptibility has a divergent singularity of the form 



with the exponent 

and the critical amplitude 



X^K{1-U/Ue)-^'^ (1.3) 



< = I (1-4) 



a; = -0.0568 ± 0.0008 (1.5) 



Several years ago, we reported some results on complex-temperature singularities of the 
susceptibility and correlation length for the Ising model @|, including a proof that for the 
square lattice, the (zero- field) susceptibility can have at most finite non-analyticities on 
the border of the complex-temperature extension of the symmetric phase, apart from its 
divergence at the physical critical point. Before proceeding, we mention some related work. 
Ref. contains some discussion of complex-temperature singularities of the Ising model on 
various 3D lattices. In a calculation and analysis of low-temperature series for the spin 1/2 
and spin 1 Ising models on the square lattice, Guttmann noted that for the spin 1/2 model x 
has a complex-temperature singularity consistent with the form (l-|-'u)~'^=, with 7^ = 3/2 (and 
u = e~'^^) 0. This paper also reported results for complex-temperature singularities for the 
spin 1 Ising model on the square lattice. Slightly later, the present authors (unaware of Ref. 
H]) carried out a similar study of the low-temperature series, inferred the same singularity 



in X at u = —1 and value of the exponent 7^ and calculated the critical amplitude [|T0| 



One of the results of Ref. |]T0| of greatest potential value for constructing conjectures for 
the exact susceptibility was a unified treatment (c.f. eqs. (3.3.5)-(3.3.9)) of both the critical 
singularity at m = 3 — 2-\/2 and the singularity at m = — 1 which yielded the relation 

7: = 2(7'-l) (1.6) 

where 7' = 7/4 is the exponent at the physical critical point, as approached from the FM 
phase. Ref. |jlO| also showed, using exact results, that universahty, in the sense that critical 



exponents are independent of lattice type, is violated at complex-temperature singular points 
and provided an explanation of this as being related to the fact that the Hamiltonian is not, 
in general, a real number at such points. This work also showed that scaling and hyperscaling 
relations are not, in general, satisfied at complex-temperature singular points. In Ref. [|lOl we 
emphasized the fact that, in accordance with the theorem of Ref. [H, X is finite (although 
non-analytic) as one approaches the point u = —1 from within the complex-temperature 
extension of the symmetric, paramagnetic (PM) phase, whereas it is divergent when one 



approaches this point from within the complex-temperature extension of the ferromagnetic 
(FM) or antiferromagnetic (AFM) phases, which behaviour is unprecedented for physical 



critical points. Hereafter, we shall denote Ref. [|10[ as I. In the present paper, we shall 
continue the study of complex-temperature singular points of the Ising model, dealing with 
the triangular lattice. A companion paper treats the closely related honeycomb lattice [|TT |. 



2 Complex Extensions of Physical Phases 

In this section we define the model and our notation, which follows that in paper I. We 
consider the Ising model on the triangular lattice (coordination number g = 6) at a tempera- 
ture T and external magnetic field H defined in standard notation by the partition function 
Z = J2{cri} e^'^^ with the Hamiltonian 

n = -JY.^^^J-HY,ai (2.1) 

where ctj = ±1 are the Z2 spin variables on each site i of the lattice, /3 = (ksT)'^, J is the 
exchange constant, < ij > denote nearest-neighbor pairs, and the magnetic moment fi = 1. 
We use the standard notation K = pj, h = PH, u as defined above in eq. ( |1.2|) , and 

v = iim\iK (2.2) 

z = e-2^ (2.3) 

which are related by the bilinear conformal mappings 

1 — "^ 1 — 2; ,„ ,, 

z = —— , V = —— 2.4 

1 + V I + z 

The reduced free energy per site is / = PF = —N~'^ limAr^^oo In Z (where Ng denotes the 
number of sites on the lattice), and the zero-field susceptibility is xo = Qfj \h=o, where 
M{H) denotes the magnetisation. Henceforth, unless otherwise stated, we consider only the 
case of zero external field and drop the subscript on xo- K is convenient to define the related 
quantity x — P~^X- After the pioneering solution for f{K, h = 0) on the square lattice by 
Onsager |T^, several papers obtained solutions for f{K, h = 0) on the triangular lattice [|T3| . 



For reviews of the Ising model on the triangular lattice, see Refs. |jTj, |T5|, |T6[. After Yang's 
derivation of the spontaneous magnetisation M for the square lattice ^7^, Potts calculated 
M for the triangular lattice |jl8| . 



We denote the critical coupling separating the paramagnetic, Z2-symmetric high tem- 
perature phase from the phase with spontaneously broken Z2 symmetry and ferromagnetic 



long-range order as Kc- As was noted above, Uc = 1/3, or equivalently, Zc = l/"\/3, and 
Vc = 2 — v^ = 0.267949.... As is well known, the Ising model on the triangular lattice (with 
equal exchange couplings along the three lattice vectors) does not have any phase transition 
to a phase with antiferromagnetic long-range order. Indeed, even at T = (corresponding 
to f = — 1 and K = z = u = oo), there is no AFM long-range order |jT9|. 

Since we shall study the susceptibility as a function of complex (inverse) temperature, it 
is useful to discuss the complex-temperature extensions of the physical PM and FM phases. 
As discussed in Ref. |P], for zero external field, there is an infinite periodicity in complex 
K under certain shifts along the imaginary K axis as a consequence of the fact that the 
spin-spin interaction (TjCTj in Ti is an integer. In particular, there is an infinite repetition of 
phases as functions of complex K] this infinitely repeated set of phases is reduced to a single 
set by using the variables f , z and/or m, owing to the fact that these latter variables have 
very simple properties under complex shifts in K: 

K ^ K + nin =^ {v ^ v, z ^ z, tt — > m} (2.5) 

K ^ K +{n + -)iiT ^ {v -^1/v, z -^ -z, u ^ u} (2.6) 

where n is an integer. On a lattice with an even coordination number q, it is easily seen 
that these symmetries imply that the magnetisation and susceptibility are functions only of 
u. Complex-temperature curves where the free energy is non-analytic were first determined 
for the square lattice by Fisher ^^. For the triangular lattice, the natural boundaries of 



the free energy, which also serve as the boundaries of the complex-temperature phases, are 
determined in a similar way as the locus of points where the argument of the logarithm 
vanishes in the (reduced) free energy^ 

f = \n2 + - / -^In cosh3(2K) + sinh3(2i^)-sinh(2i^)P(^i,^2) (2.7) 



where 



P{ei, 62) = cos Oi + cos 02 + cos(^i + 62) (2.8) 

Expressed as a function of u, this is 

f = 3K + -J_^ J ^ -^ ln[(l + 3u') - 2u{l - u)P{e^, 62)] (2.9) 

The function P{9i,62) ranges from a maximum value of 3 at ^1 = 6*2 = to a minumum 
value of —3/2 at 61 = 62 = 27r/3. Hence the above curves are the locus of solutions to the 

^ The free energy is trivially infinite aX K = 00; since this is an isolated point and hence does not form 
part of a natural boundary separating phases, it will not be important here. 



equation 



3 

l + 3u^ - 2m(1 -u)x = for --<x<3 (2.10) 



2 



This locus consists of the union of the circle 



u = -l + le'^ (2.11) 



3 3 

for < (p < 27T with the semi-infinite hne segment 

1 
-oo<u< — (2.12) 

as pictured in Fig. 1(a). We label the endpoint (e) of this semi-infinite line as Ue, given 
above in eq. ( |1 . 1| ) . The intersection of the circle ( |2.11|) and the line segment ( |2.12| ) consists 



of the single point 

Us = -1 (2.13) 

which will also be important in our later discussion. Since eq. (|2.10|) has real coefficients, the 
solutions are either both real or are complex conjugate pairs, which explains the reflection 
symmetry of Fig. 1(a) about the Im{u) = axis. For x = 3, eq. (|2.1CI|) has a double root at 
the physical critical point u = Uc = 1/3. As x decreases from 3 to 0, the complex conjugate 
roots move from u = 1/3 along the circle (|2.11| ) to ±i/^/3. As x decreases from to —1, the 



roots move from these points around the circle and finally join again to form a double root 
at -u = — 1 for X = — 1. The behaviour is qualitatively different as x decreases from —1 to its 
minimum value —3/2: in this interval, the solutions consist of two real roots. These move 
outward from m = — 1 along the negative real axis. As x ^ —3/2, one root moves rightward 
to u = —1/3 while the other moves leftward to u = — oo. 

The curves in Fig. 1(a) divide the complex u plane into two regions which are complex- 
temperature extensions of physical phases: (i) the complex FM phase, which is the extension 
of the physical FM phase occupying the interval < u < Uc, and (ii) the complex PM phase, 
which is the extension of the physical PM phase occupying the interval Uc < u < 1. As 
is evident from Fig. 1(a), there is a section of the line segment ( |2.12| ) protruding into the 
interior of the complex FM phase. All directions of approach to the endpoint Ue of this line 
segment are from within the complex FM phase, except for the approach precisely from the 
left, along the line segment itself. Note that the outer part of the line segment extending 
to u = — oo prevents one from making a circuit around the origin staying within either the 
complex FM or PM phase. 

This suffices for our analysis of x, since it is only a function of u. However, especially for 
our study of complex-temperature singularities of the Ising model on the honeycomb lattice 



in a companion paper, it is also useful to exhibit the corresponding curves in the variable z 
and V. The former are directly obtained from eqs. (|2.11| )- (|2.12| ) and are shown in Fig. 1(b). 
One observes that, whereas the complex w-plane consists of two phases, the complex z-plane 
is divided into three phases; in addition to the complex-temperature extensions of the PM 
and FM phase, there is also a third, denoted O (for "other") in Fig. 1(b) lying outside of the 
boundary of the complex FM phase and to the left of the semi-infinite vertical lines running 
upward from Ze+ = z/v^ to zcxd and downward from z^,- = —i/\/3 to —ioo. The physical 
phases are, of course, the (i) FM phase, in the interval < z < Zc = l/y/S and (ii) the PM 
phase, in the interval Zc < z < oo. The O phase has no overlap with any physical phase. 
The same type of phenomenon was true of the complex-temperature phases for the square 
lattice (see Fig. 1(b) in our previous paper I). 

In order to obtain the analogous boundary curves in the v plane, one may either use the 
bilinear conformal mapping ( p.4|) on the curves in Fig. 1(b) or express the free energy in 



terms of v and solve for the points where the argument of the logarithm in the integrand 
vanishes. One finds that, in addition to the isolated point v = 1 {K = oo), / is non-analytic 
along the locus of points satisfying the equation 

l-2t; + 6t;2-2t;^ + t;^-2t;(l-t;)2a; = - 3/2 < a; < 3 (2.14) 

where, as above, x = P{9i,92). Since this is a quartic equation with real coefficients, the 
solution for each x consists generically of two pairs of complex-conjugate roots (which may 
include real double roots). Furthermore, under the mapping v —>■ 1/v, the left-hand side 
of eq. (|2.14|) transforms into itself, up to an overall factor of f ~^. This implies that the 



solutions of this equation occur in reciprocal pairs. We can describe the solutions as follows. 
For X = 3, there are two solutions: the physical critical point f c = 3 — 2v^, and its inverse, 
1/vc = 3 + 2v2. As X decreases through the range from 3 to —1, the double root at f = 1/vc 
splits into a conjugate pair, the members of which move along the outer curves in Fig. 1(c) 
from 1/vc over to i and —i, respectively, while the double root at Vc splits into another 
complex conjugate pair, the members of which move along the inner curve in this figure, to 
—i and i. For x = —1, two roots coalesce at i and the other two at —i to form two double 
roots. As X decreases from —1 to its minimum, —3/2, these double roots split into four roots, 
which move outward from ±i along the unit circle \v\ = 1. At x = —3/2, the rightward- 
moving roots reach the values exp(±z7r/3), while the leftward- moving roots coalesce into a 
single double root at f = —1. The physical phases are (i) PM, for < v < Vc, and (ii) 
FM, for fc < f < 1. The regions marked PM and FM in Fig. 1(c) are the extensions of 
these to complex-temperature. The outermost region in the figure, marked O is a phase 
which has no overlap with any physical phase. The point u = Us = —1, or equivalently, 

6 



the points z = ±i, v = ^i are singular points of the curves, in the terminology of algebraic 
geometry ||21||; specifically, they are multiple points of index 2. The singular points of the 
natural boundaries for the Ising model on the square lattice also occur at these points and are 
also of index 2. However, the topology of the complex phases is different for the triangular 
and square lattices, reflecting the fact that the former is a tight-packed lattice without an 
AFM phase, while the latter is a loose-packed lattice. 

As noted in paper I, using the general fact that the high-temperature and (for discrete 
spin models such as the Ising model) the low-temperature expansions have finite radii of 
convergence, we can use standard analytic continuation arguments to establish that not just 
the free energy and its derivatives, but also the magnetisation and susceptibility, are analytic 
functions within each of the complex-temperature extensions of the physical phases. This 
defines these functions as analytic functions of the respective complex variable {K, v, z, u, 
or others obtained from these). We refer the reader to paper I for other relevant background, 
definitions, and notation. 

3 Analysis of Low- Temperature Series 

3.1 Generalities 

The low-temperature series expansion for x on the triangular lattice is given by 

oo 

X = 4^3(1 + ^c„m") (3.1.1) 

ra=l 

This expansion has a finite radius of convergence and, by analytic continuation from the 
physical low-temperature interval < u < Uc, applies throughout the complex extension of 
the FM phase. Since the u^ factor is known exactly, it is convenient to study the remaining 
factor Xr = u~^x- For the triangular lattice, following earlier work P], the expansion coef- 
ficients Cn,t were calculated by the King's College group to order n = 13 (i.e., x to 0{u^^)) 
in 1971 P2[, and to order n = 18 (i.e., x to 0(m^^)) in 1975 [El]. As discussed above, it 



was noticed quite early that the low-temperature series does not give very good results for 
the position of the physical critical point or for its exponent (as this point is approached 
from within the FM phase), 7'. The cause for this was recognized to be the existence of 
the unphysical, complex-temperature singularity at Me, the same distance from the origin 
as the physical singularity at tic H- We see that this is precisely the endpoint of the 
semi-infinite line segment discussed in section ^ which protrudes into the interior of the 
complex FM phase. Since at that time a primary goal was to study the physical critical 



point and, e.g., test the lattice universality of critical exponents and various scaling rela- 
tions, the complex-temperature singularity was understandably regarded as a nuisance, and 
methods were devised to transform the series expansion variable so as to map it farther from 
the origin than the physical critical point. We have noted above that Guttmann, continuing 
his earlier work on complex-temperature singularities of the Ising model on various lattices 
0, H, carried out a study of this singularity at m = —1/3 on the triangular lattice and 
determined the corresponding exponent and critical amplitude, using ratio and dlog Pade 
methods. 

Here we shall re-analyse this singularity using, in addition to dlog Pade approximants the 
powerful method of differential approximants. (For reviews of these methods, see Refs.[^], 
p^ , and p5|.) Given a knowledge of the locus of points where the free energy is non- analytic 



and where the magnetisation vanishes continuously, we have also been motivated to explore 
the behaviour of x ^^ the point m^ = — 1. Of course, a first thing to check is whether the 
low-temperature series for x o^i the triangular lattice has been calculated to higher order 
since the King's College work in 1975; we have found that apparently it has not p6|, ^ . 
As one approaches a generic complex singular point denoted sing from within the complex- 
temperature extension of the FM phase, x is assumed to have the leading singularity 

X{U) ~ A^„^(l - u/UsingY^''"''^^ (l + ai(l - u/Using) + •••) (3.1.2) 

where A'^^^g and 7^^^^ denote, respectively, the critical amplitude and the corresponding 
critical exponent, and the ... represent analytic confluent corrections. One may observe that 



we have not included non-analytic confluent corrections to the scaling form in eq. (|1.3|) . 
The reason is that, although such terms are generally present at critical points in statistical 
mechanical models, previous studies have indicated that they are very weak or absent for 
the usual critical point of the 2D Ising model p^, ^. The dlog Pade and differential 
approximants then directly give Using and "^'sing- As noted above, since the pref actor v? is 
known and is analytic, we actually study Xr- 

3.2 Exponent 7^ at w = —1/3, Revisited 

We first re-examine the singularity at u = —1/3 previously studied in detail by Guttmann 
||T|]. This singularity is studied by an analysis of the low-temperature series expansion for x, 
since the approach to this point from all directions in the complex u plane except directly 
from the left along the singular line segment ( |2.12| ) are from with the (complex-temperature 
extension of the) FM phase. For exploratory purposes, we first carried out a simple dlog 
Pade study of the series for Xr in the variable u. While this yielded a reasonably stable 



location for the singularity, it did not yield very stable or accurate values for the exponent. 
An obvious reason for this is the presence of the physical critical point the same distance 
away from the origin. Accordingly, we transformed the series to one in the variable 

u 

which has the effect of mapping the physical critical point to u = cxd and the complex- 
temperature singularity at Ue to Ue = —1/6 = —0.16666.... A dlog Pade study of this series 
gives improved results for the location and exponent of the singularity. We obtained by far the 
most precise results using differential approximants. As is well known, this method is more 
powerful than dlog Pade approximants since it is able to represent an analytic background 
term as well as a singular term in a given function. In this method, the function / = Xr(C) 
being approximated satisfies a linear ordinary differential equation (ODE) of K'th order, 
CM,LfK{0 = Ef=oQAOD'fK{0 = R{0, where Q,{0 = eSo Q,X' and R{C) = EtoRiC 
(and ( denotes a generic complex variable). We shall use the implementation of the method 
in which D = Qd/dC,. The solution to this ODE, with the initial condition /(O) = 1, is the 
resultant approximant, labelled as [L/Mq; ...; Mk\- The general solution of the ODE has the 
form /j(C) ~ ^i(C)IC ~ Cjl"^'' + -S(C) foi^ C ~^ Cj- The singular points Q are determined 
as the zeroes of Qk{C) aiid are regular singular points of the ODE, and the exponents are 
given by —pj = K — 1 — Q k-i{Cj) / iCjQ'KiCj)) ■ Given the above-mentioned evidence that for 
the 2D Ising model, non-analytic confluent singularities are weak or absent for the physical 
critical point, which motivates the form ( |1.3D , K = 1 differential approximants should be 
sufficient for our purposes here. (Of course, K = 2 differential approximants could be used 
in a further study to investigate confluent singularities.) We also choose to use unbiased 
differential approximants since we wish to confirm that the location of the singularity is, 
indeed, at m = —1/3 as expected. The results of this analysis are given in Table 1. As 
expected, we obtain higher precision results using differential approximants on the series for 
Xr in the transformed variable u. These are given in Table 2. 

The results of this study agree with the old inference of a singularity in % at m = Wg = 
— 1/3 (i.e., u = —1/6). We have found that the differential approximants yield a consid- 
erably more precise determination of the exponent than the dlog Pade approximants, and 
our results are in excellent agreement with Guttmann's earlier conclusion using ratio and 
dlog Pade methods that j'^ = 5/4 |]T|. Specifically, we find 

7^ = 1.249 ± 0.005 (3.2.2) 

As will be discussed in our companion paper on the honeycomb lattice 0], the result (|3.2.2|) 



is closely related to two findings which we report there, viz., that for the Ising model on the 



[L/Mo; Ml] 


Using 


I'^sing ^e|/ "^el 


< 


[0/7 


7] 


-0.3332767 


1.7 X 10-^ 


1.2426 


0/7 


8] 


-0.3333619 


0.86 X 10-^ 


1.2522 


0/7 


9] 


-0.3333237 


2.9 X 10^5 


1.2476 


0/8 


7] 


-0.3333691 


1.1 X 10-^ 


1.2530 


0/8 


8] 


-0.3333447 


3.4 X 10"^ 


1.2502 


0/9 


7] 


-0.3333107 


0.68 X 10-^ 


1.2460 


1/6 


6] 


-0.3334340 


3.0 X 10-^ 


1.2594 


[1/6 


7] 


-0.3332574 


2.3 X 10-^ 


1.2407 


[1/6 


8] 


-0.3333707 


1.1 X 10-^ 


1.2534 


[1/7 


6] 


-0.3332608 


2.2 X 10-^ 


1.2410 


[1/7 


7] 


-0.3333296 


1.1 X 10-5 


1.2478 


[1/7 


8] 


-0.3334177 


2.5 X 10-^ 


1.2602 


[1/8 


6] 


-0.3333782 


1.3 X 10-^ 


1.2542 


[2/6 


6] 


-0.3334040 


2.1 X 10-^ 


1.2581 


[2/6 


7] 


-0.3333311 


0.68 X 10-5 


1.2494 


2/6 


8] 


-0.3333613 


0.84 X 10-^ 


1.2524 


2/7 


6] 


-0.3333310 


0.69 X 10-5 


1.2494 


2/7 


7] 


-0.3333301 


0.96 X 10-5 


1.2484 


2/8 


6] 


-0.3333618 


0.85 X 10-^ 


1.2525 


3/4 


6] 


-0.3333208 


3.8 X 10-5 


1.2282 


[3/6 


6] 


-0.3332486 


2.5 X 10"^ 


1.2374 


[4/6 


4] 


-0.3333540 


0.62 X 10-^ 


1.2553 


[5/5 


6] 


-0.3332986 


1.0 X 10"^ 


1.2423 


[6/4 


6] 


-0.3333335 


0.47 X 10-6 


1.2463 


[6/5 


5] 


-0.3332566 


2.3 X 10-^ 


1.2375 


[6/6 


4] 


-0.3334199 


2.6 X 10-^ 


1.2624 


[7/4 


4] 


-0.3333719 


1.2 X 10-^ 


1.2570 


[8/4 


3] 


-0.3333846 


1.5 X 10-^ 


1.2596 


8/4 


4] 


-0.3333423 


2.7 X 10-5 


1.2511 



smg and 7e from differential approximants to low-temperature series for 



Table 1: Values of u 

Xr{u). See text for definition of [L/Mq,Mi 

satisfy the accuracy criterion \usinn — Up\/\u 



I'^sing "•€ 



approximant. We only display entries which 
< 3 X 10"^ where u^ = -1/3. 
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[L/Mo; Ml] 


"^sing 


\'^sing "^el/ "^el 


ie 


[0/6 


7] 


-0.1666760 


0.56 X 10-^ 


1.2516 


[0/6 


8] 


-0.1666666 


1.0 X IQ-^ 


1.2481 


0/7 


6] 


-0.1666760 


0.56 X 10-^ 


1.2517 


0/7 


7] 


-0.1666643 


1.4 X 10-^ 


1.2471 


0/7 


8] 


-0.1666621 


2.7 X 10-5 


1.2462 


0/7 


9] 


-0.1666628 


2.3 X 10-5 


1.2466 


[0/8 


6] 


-0.1666666 


4.9 X 10"^ 


1.2481 


[0/8 


7] 


-0.1666621 


2.7 X 10-5 


1.2462 


[0/8 


8] 


-0.1666629 


2.3 X 10-5 


1.2466 


[0/9 


7] 


-0.1666628 


2.3 X 10-5 


1.2466 


[1/6 


6] 


-0.1666543 


0.74 X 10-^ 


1.2424 


[1/6 


7] 


-0.1666587 


0.48 X 10-^ 


1.2443 


[1/6 


8] 


-0.1666584 


0.50 X 10-^ 


1.2441 


1/7 


6] 


-0.1666587 


4.8 X 10-5 


1.2443 


1/7 


7] 


-0.1666582 


5.1 X 10-5 


1.2440 


[1/7 


8] 


-0.1666749 


4.9 X 10-5 


1.2535 


[1/8 


6] 


-0.1666583 


5.0 X 10-5 


1.2441 


1/8 


7] 


-0.1666748 


4.9 X 10-5 


1.2535 


2/5 


7] 


-0.1666817 


0.90 X 10-^ 


1.2564 


[2/6 


6] 


-0.1666742 


0.45 X 10-^ 


1.2524 


[2/6 


7] 


-0.1666577 


0.54 X 10-^ 


1.2438 


[2/6 


8] 


-0.1666630 


2.2 X 10-5 


1.2467 


[2/7 


5] 


-0.1666817 


0.90 X 10-^ 


1.2565 


[2/7 


6] 


-0.1666577 


0.54 X 10-^ 


1.2438 


[2/7 


7] 


-0.1666619 


2.8 X 10-5 


1.2461 


[2/8 


6] 


-0.1666631 


2.2 X 10-5 


1.2467 



Table 2: Values of Using and 7^ from differential approximants to low-temperature series for 
Xr{u). See text for definition of [L/Mo,Mi] approximant. We only display entries which 
satisfy the accuracy criterion 



u 



sing 



Up 



\U. 



< 1 X lo- 



ll 



[3/5 


6] 


-0.1666780 


0.68 X 10"^ 


1.2545 


[3/5 


7] 


-0.1666593 


0.44 X 10-^ 


1.2446 


[3/6 


5] 


-0.1666780 


0.68 X 10"^ 


1.2545 


[3/6 


6] 


-0.1666564 


0.61 X 10-^ 


1.2431 


[3/6 


7] 


-0.1666615 


3.1 X 10-5 


1.2458 


3/7 


5 


-0.1666593 


0.44 X 10"^ 


1.2446 


3/7 


6 


-0.1666615 


3.1 X 10-5 


1.2458 


4/5 


5 


-0.1666748 


0.49 X 10-^ 


1.2528 


4/5 


6 


-0.1666575 


0.55 X 10-4 


1.2436 


4/5 


7 


-0.1666632 


2.1 X 10-5 


1.2468 


4/6 


5 


-0.1666575 


0.55 X 10-4 


1.2436 


[4/6 


6] 


-0.1666593 


0.44 X 10-4 


1.2444 


[4/7 


5] 


-0.1666633 


2.0 X 10-5 


1.2469 


[5/5 


5] 


-0.1666577 


0.54 X 10-4 


1.2437 


[5/5 


6] 


-0.1666634 


1.9 X 10-5 


1.2470 


[5/6 


4] 


-0.1666809 


0.85 X 10-4 


1.2573 


[5/6 


5] 


-0.1666637 


1.8 X 10-5 


1.2472 


[6/4 


6] 


-0.1666618 


2.9 X 10-5 


1.2458 


[6/5 


3] 


-0.1666838 


1.0 X 10-4 


1.2594 


6/5 


5 


-0.1666630 


2.2 X 10-5 


1.2467 


[6/6 


4 


-0.1666672 


3.2 X 10-6 


1.2496 


[7/4 


5 


-0.1666595 


0.43 X 10-4 


1.2444 


7/5 


4 


-0.1666691 


1.5 X 10-5 


1.2507 


[8/5 


3 


-0.1666811 


0.87 X 10-4 


1.2578 



Table 3: continuation of Table 2. 
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honeycomb lattice, the (uniform) susceptibihty x and the staggered susceptibihty x*-"^ have 
divergent singularities at the point z = zii = —1 with respective exponents which we infer to 
have the values 7^ = 7^ a = 27^ = 5/2. 

3.3 Critical Amplitude at w = —1/3 Singularity 

Here we have used a method complementary to that employed earlier by Guttmann to extract 
the critical amplitude A'^ [|T|. Specifically, we compute the series for {x^)^^'^''■ Since the exact 
function (Xr)^^^" has a simple pole at Ug, one performs a Pade analysis on the series itself 
instead of its logarithmic derivative. The residue at this pole is 

R, = -u,{A',^,f'^'^ (3.3.1) 

where A'^^ denotes the critical amplitude for Xr- It follows that 

K = ^<K,e = -m~"'Rl" (3.3.2) 

Using the inferred value 7g = 5/4 to compute the series, extracting the residue i?e, and then 
using eq. ( p.3.2| ), we obtain the critical amplitude. We have also done the same analysis on 



the transformed series in the variable u = u/{l — 3m) for which the physical critical point is 
mapped to infinity, and Ue to Ue = —1/6. Equating the singular forms 

X{u) ~ A;(1 - n/«,)-^/^ 

= x{u)r^K{l-u/ue)~'^^ (3.3.3) 

we have 

A'^ = 25/^i'g (3.3.4) 

In our analysis of the Pade results, we plot the individual determinations of the critical 
amplitude from each [L/M] approximant as a function of the accuracy with which that 
approximant locates the singular point, as measured by the relative deviation \using — Ue\/\ue\ 
and \using — Ue\/\ue\- Wc obtain the final estimate of the critical amplitude by extrapolating 
these results to zero deviations. Combining our results from both analyses, we obtain the 
result 

a; = -0.05766 ± 0.00015 (3.3.5) 

This agrees with the result by Guttmann [|l|, cited in eq. (|1.5| ) and has a somewhat smaller 
estimated uncertainty. 
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It is of interest to compare this with the critical amphtude at the physical critical point, 
as approached from within the FM phase. For reference, we note that many authors express 
the singularity in terms of T rather than u, viz., 

~ A'^il-u/u,)-''' (3.3.6) 

with 7' = 7/4. The critical amplitudes are related according to 

A', = m,y/'K,T (3.3.7) 

An early analysis by Essam and Fisher of a low-temperature series for x up to 0(m^°) for x 
yielded the resuh [||] A^^ = 0.0248 ±0.0006, or equivalently, A'^ = 0.0292 ±0.0007. A more 
accurate determination was made by Guttmann 133], based on the series for x calculated 



to order ©(m^^) H, with the result A[^ = 0.0246 ± 0.0002, i.e., A^ = 0.0290 ± 0.0002. In 



the same paper, using a generalised law of corresponding states, as formulated in Ref. |pO | 



(further discussed in Ref. [0), Guttmann obtained a high-precision analytic determination 



of the critical amplitudes for the Ising susceptibility for the triangular, honeycomb, and 



kagome lattices (analogous to that obtained in Ref. [28| for the square lattice). For the 
triangular lattice, he calculated 

A'^T = 0.0245189020 i.e., A'^ = 0.028905388 (3.3.8) 

Using this in conjunction with eq. ( p.3.5|) , we find 



A' 

-f = -(1.995 ±0.005) (3.3.9) 

where the uncertainty arises completely from the uncertainty in A'^. One may observe that 
the ratio (p.3.9|) is slightly less than, but close to, the simple relation A'^/A'^ = —2. We should 
note, however, that we do not know of a specific analytic basis for such a simple amplitude 
relation. 

3.4 On the Singularity at u = —1 

Although the complex-temperature phase diagram of the Ising model on the triangular lattice 
is different from that on the square lattice, the point u = Ug = —1 is somewhat similar for 
both of these lattices, in the sense that at this point on the two respective lattices, (i) 
the spontaneous magnetisation and the inverse correlation length vanish continuously as 
one approaches this point from within the complex FM phase (as will be discussed further 
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below). The differences include the fact that (ii) on the square lattice, the complex PM, 
FM, and AFM phases are all contiguous at m = — 1, whereas on the triangular lattice, only 
the complex FM and PM phases are contiguous (there is no AFM phase); and (iii) on the 
square lattice, the point u = —1 can be reached by moving in a straight line outward along 
the negative real axis from the origin in the w-plane, whereas on the triangular lattice, it 
cannot be; one first encounters the endpoint of the singular semi-infinite line segment ( p.l2|) 
at u = —1/3. It is thus naturally of interest to investigate the behaviour of the susceptibility 
at M = —1. 

For this purpose, we have carried out an analysis of the low-temperature series for x 
on the triangular lattice, using dlog Pade and differential approximant methods. However, 
we have found that the series, at least at the order to which it has been calculated, cannot 
probe this point very sensitively. We believe that the reason for this is the property (iii) listed 
above, viz., the fact that the series is dominated by the singularity at m = —1/3, which is not 
only closer to the origin but also directly in front of the point u = — 1 as approached from 
the origin. This is just as true of the series transformed to the variable u as it is of the series 
in u since the image of the point u = —1/3 is at -S = —1/6, directly in front of the image of 
the point u = —1, at -u = —1/4, as approached from the origin in the u plane. Hence, one 
does not expect the series transformed to the u variable to be more sensitive to the u = —1 
singularity, and, indeed, it is not. We have analysed the series with many transformations 
of expansion variables which are designed to map both the u = —1/3 singularity and the 
physical critical point at m = 1/3 farther away from the origin than the u = —1 singularity. 
Examples include Utrans. = sinh('u)/(l — 9u^) and Un-ans. = Mexp(a(l + -u)^ — a) for various 
a > and positive even integer p. However, even in the transformed variables, the series 
yield inconclusive results for 7^. We can say that if the scaling relation a' + 2/3 + 7' = 2 holds 
at M = — 1 for the Ising model on the triangular lattice, as it does on the square lattice, then, 
given the exact result that Ps = 3/8 (see eq. ([4.1.1|) below) and our finding that a'^ = from 
the exact free energy (see eq. ( |6.3.1|) below), it would follow that 7^ = 5/4, which would be 
equal to the value of the exponent 7^ for the singularity at Ue- 

In passing, we note that the analogue of the transformation to an elliptic modulus variable 
which we found to be very valuable in the case of the square lattice is not helpful here, for 
reasons which can easily be explained. The requisite elliptic modulus is clear from the exact 
solutions for the free energy []T3[ and spontaneous magnetisation \TB\ of the model (the latter 
of which can be written as M = (1 — (A;<)^)^/^) and is given by 

<" (1 + 3m)i/2(i_^)3/2 (3.4.1) 

First, the image of the singular boundary ( p.llj )-( p.l2| ) extends all the way in to the origin 
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in the complex /c<-plane. Specifically, the semi-infinite line segment ( p.l2|) is mapped via a 



2- valued homomorphism to the segment — 1 < (k^) < 0. This is quite different from the 
situation for the square lattice, where the mapping from u to the elliptic modulus variable 
for that lattice, viz., k^ = 4^/(1 — u)"^, mapped the boundary of the complex FM phase to 
the unit circle |/c<| = 1, so that no portion of this boundary impinged upon the origin in the 
complex (fc<)-plane. A second important difference is that whereas for the square lattice, 
|fc<| is less than unity in the complex FM phase, for the triangular lattice, |/c<| is not only 
not less than unity in this phase; it actually diverges as one approaches the endpoint Ue of 
the singular line segment protruding into this phase. To see this, let 



u = -l + \ee''^ (3.4.2) 



where e is real and positive. Then, 



k< -^ ; — . ,,-, ,„ as e — >• (3.4.3) 

^ 2(ee*<^)V2 ^ > 



Taking the branch cut for the fractional powers in (|3.4.1| ) to lie from u = Uetou = — oo, then 



for the approach iou = —1/3 from the origin along the negative real axis, which corresponds 
to = in eqs. (|3X2D and (|3A^ ), k< {i/2){e)-^/^ -^ -too. 



4 Behaviour of Magnetisation and Implications 

4.1 Divergence of M at w = —1/3 

In discussing our results for the singularity in % at m = -Ue — 1/3, we first note that some 
insight into the location of this singularity can be gained from an examination of the exact 
expression for the spontaneous magnetisation of the Ising model on the triangular lattice 

M = c-±^y" C-^y" (4.1.1) 

\l-uj V1 + 3m/ ^ ^ 

which applies within the physical FM phase and, by analytic continuation, throughout the 
complex-temperature extension of this phase as shown in Fig. 1(a). (Of course, outside of 
this complex-extended FM phase, M vanishes identically.) M diverges^ at the point Ue, with 
exponent 

Pe = -l (4.1.2) 



^ For physical temperatures, the unit-normahzed magnetisation is bounded to he between and 1. 
However, this bound does not apply for complex temperatures, and, indeed, nothing prevents M from 
having divergent singularities for complex temperatures. 
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Thus it is not surprising that at a point such as Ue where there is a singularity in M, there 
is also a related singularity in x- 

4.2 Theorem onM^oo^X^cxD 

Such an exotic phenomenon as a divergent spontaneous magnetisation has received very little 
attention in the literature. We begin by stating and proving a theorem which deals with an 
important effect of such a divergence. First, we prove a lemma concerning 2-spin correlation 
functions: 

Lemma 1 

Assume that a given statistical mechanical model has a phase with ferromagnetic long-range 
order. In this phase, and in its extension to complex temperature, the 2-spin correlation 
function can always be written in the form 

< cr^an' >= M^c{n,n') (4.2.1) 

where M is the spontaneous magnetisation and c{n, n') contains all of the dependence on 
the lattice sites n and n'. 

Proof 

This result follows easily from the fact that one of the equivalent definitions of M is precisely 
via the relation 

M^ = lim < a„cr„/ > (4.2.2) 

Given the correlation function < cr„cr„/ >, one may thus calculate M^ from the limit (|4.2.2|) 
and divide, thereby obtaining the function c{n,n') (with the property lim|„„„/|_^oo c(n,n') = 
1). D 
As an immediate corollary, we have 

Lemma 2 

Assume that a given statistical mechanical model has a phase with ferromagnetic long-range 
order. In this phase, and in its extension to complex temperature, the connected 2-spin 
correlation function can always be written in the form 

< a„Cr„/ >conn.= M"^ c{n,n')conn. (4.2.3) 

where M is the spontaneous magnetisation and c(n, n')conn. contains all of the dependence 
on the lattice sites n and n' . 
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Proof 

This follows immediately from the definition of the connected 2-spin correlation function as 
< cTnan' >conn. = < cr„cr„/ > — M^, which also shows that c(n, n')conn = c{n, n') — 1. □ 
We then proceed to 

Theorem 1. 

If the magnetisation M diverges as one approaches a given point from within the (complex- 
temperature extension of the) ferromagnetic phase, then, in the same limit, the susceptibility 
also diverges. 

Proof 

The susceptibility x is given as the sum over the connected 2-spin correlation functions 

X = N~^ ,lim J2 < ^«^«' >conn. 
n,n' 

= ^ < (ToCr„ >conn. (4.2.4) 

n 

(where the homogeneity of the lattice has been used in the second line). Using Lemma 2, 
we have 

X = M'J2<ri,n%onn. (4.2.5) 

n 

It follows that, in general, a divergence in M will cause a divergence in x. □ 
Note that this would be true even in the hypothetical case in which M is divergent but the 
correlation length is finite, so that the sum J2n c{n, n')conn. converges. In passing, we note 
that if M diverges at a given point uq (or the corresponding points zq), this will correspond 
to the infinite set of values of K which map to uq or the set of zq according to the symmetries 
( p.5| ) and, for even q, ( p.6| ). 



We next apply Theorem 1 to our immediate topic of study: 
Corollary 1. 

For the Ising model on the triangular lattice, x has a divergent singularity at u = —1/3. 
Proof 



This follows from the fact that, as we know from the exact result, (^4.1.1|) , M (analytically 
continued throughout the complex-temperature extension of the FM phase) diverges at u = 
—1/3 together with Theorem 1. □ 

Note that unlike the study of the low-temperature series, which is, of course, approximate 
since the series only extends to finite order, this is an exact rigorous result. What the studies 
by Guttmann [|l| and the present authors yielded beyond the result of the theorem was the 
actual values of the exponent 7g and critical amplitude A'^. 

Also, observe that we did not explicitly use any property of the correlation length to 
make this conclusion. Our theorem and corollary also allow us to infer without a direct 
calculation that the correlation length does in fact diverge at u = —1/3 (as this point is 
approached from the interior of the complex-temperature extension of the FM phase, i.e. 
from all directions except from the left along the singular line segment (|2.12| )). We can 



deduce this because if the correlation length were finite, then the sum over 2-spin correlation 



functions in (|4.2.4| ) would be finite. (Since only the large-distance behaviour is relevant to a 
possible divergence, one can replace the sum by an integral, and because of the exponential 
damping from < aQOr >~ r~^e'"/'», the integral is finite.) But then since the only divergence 
would arise from the prefactor of M^, we would have the exponent relation 7g = —2(3e- Since 
— 2/5e = 1/4, while the series analyses yield 7^ = 5/4, the above exponent relation does not 
hold. This shows then, that the correlation length diverges at m = —1/3. 

4.3 Other Singularities of M 

In addition to its divergence at u = —1/3, M as given by ( |4.1.1D vanishes continuously at 
the complex-temperature point u = Ug = —1, with exponent 

Ps = \ (4.3.1) 

as well as the usual critical point, u = Uc = 1/3 (with exponent j3 = 1/8). With the 
exception of these two points, Uc and Us, M vanishes discontinuously elsewhere along the 
outer boundary of the complex-temperature extension of the FM phase. Note that despite 
the (1 — m)~^/^ factor in the exact expression (^4.1.1|) , M does not in fact diverge at m = 1, 
since this point lies outside the region where eq. ([4.1.1| ) applies (which, we recall, is the 
region which can be reached by analytic continuation from the physical FM phase, i.e. the 
complex-temperature extension of the FM phase). 
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5 Some Remarks on Scaling Relations 

We have previously shown that in general at complex temperature singularities, a number of 
the usual scaling relations applicable for physical critical points do not hold. In particular, 



for the Ising model on the square lattice we have proved g, ^ that as one approaches the 
boundary of the complex-temperature extension of the symmetric, PM phase, the suscepti- 
bility has only finite non-analyticities, except for its divergence at the physical critical point. 
Our proof was based on the fact that in this limit, the inverse correlation length only vanishes 
at the physical critical point. From our theorem, it follows that at the points w = ±z (i.e., 
the point m = m^ = — 1), as approached from within the complex PM phase, the exponent 
describing the finite non-analyticity in x is negative: 7^ < 0. On the other hand, at this 
point Us = —1, as approached from within the interior of the complex-temperature extension 
of the FM phase, x does diverge P, |10|- As we noted in Ref. |T^, this implies, among other 



things, that at this singularity, u = m^ = — 1, 7s 7^ 7s- This violates the exponent scaling re- 
lation (derived from the usual renormalization-group method) applicable at physical critical 
points, that 7' = 7. 



Secondly, we have shown |T0[ by the use of exact results for the magnetisation, that 
universality, in the sense that critical exponents are independent of lattice type, is, in general, 
violated at complex-temperature singularities. 

Thirdly, we have shown [|l^] that for the square lattice, as one approaches the point 



u = Us = —1 from within the complex FM (or AFM) phase, the inverse correlation length 
^FM row defined from the row (or equivalently, column) connected 2-spin correlation functions 
(analytically continued throughout the complex-temperature extension of the FM phase) 
vanishes like |1 -|-m|''= with u'^ = 1, whereas the inverse correlation length ip]^^ defined from 
the diagonal connected 2-spin correlation functions vanishes like |l+ti|~''s,dias. -v^^jth Vs^iag. = 2 
in the same limit. As noted in paper I, this situation is unprecedented for physical critical 
points. Indeed, more generally, using the exact calculation of the asymptotic form of the 
2-spin correlation function < (To,oo^m,n > @] for the square lattice (analytically continued 



throughout the complex FM phase), we have shown that the inverse correlation length 
extracted from < cro,oO"m,n >conn. as r = [m? + ra^)^/^ -^ 00 vanishes as m — > m^ from within 
the complex FM phase like \u + Ij''* with z/^ = 1 if 6' = arctan(?7i/n) does not represent a 
diagonal of the lattice, i.e. is not equal to ±7r/2 or ±37r/2. This result means that for the 
Ising model on the square lattice, as one approaches u = —1 from within the complex FM 
phase, the spin-spin correlations along the diagonal become long-ranged faster than those in 
other directions. This, in turn, undermines the naive use of renormalization group methods 
to derive scaling relations for exponents since these methods rely on the existence of a single 



20 



diverging length scale provided by the correlation length. 

Nevertheless, one does find f^, |l^ that for the square lattice, as one approaches the 
point u = —1 from within the complex FM phase, the exactly known exponents a'g = 
(log divergence in specific heat) and Ps = 1/4, together with the inferred exact result @, |10 



7^ = 3/2, satisfy the scaling relation a'^ + 2/3^ + 7^ = 2. However, as we noted IT^, the 



corresponding relation for the approach to m = — 1 from within the symmetric phase, viz., 
as + 2Ps + 7s = 2 is false, since a^ = (log divergence) and, as mentioned above, 7^ < 0. 

These findings show that universality, scaling, and exponent relations which were appli- 
cable to physical critical points do not, in general, hold for complex-temperature singular 
points. In particular, for the subset of these complex-temperature singularities where M 
diverges, there is yet another unusual implication. Let us define 

t=\T-T,\/T, (5.1) 

As t — *> 0, the connected 2-spin correlation function approaches the critical decay 

< CToar >conn.^ r-'^d-2+,) (5.2) 

In the FM phase, this correlation function can also be written in the form ( [4.2.3|) , i.e., 
< o"„cr„/ >conn.= M'^c{n,n'). Then according to standard assumptions of scaling and the 
renormalization group, as t — > and r ^ oo with 

P = r/i (5.3) 

fixed, where ^ ~ t~" , the spatially-dependent term c(n,n') becomes only a function of the 
scaled length, p, i.e., 

c(n,n')~p-^''-'+^^ (5.4) 

so that in this limit, r-('^-2+'') ~ t'^'^ / (rf' Y'^-'^+'^\ and hence 

2p = u'{d-2 + ri) (5.5) 

which is one of the well-known scaling relations among critical exponents. Now recall that 
this use of the renormalization group requires a divergent correlation length, so u' must 
be positive.0 But ii (3 < 0, as is the case for a complex-temperature singularity where M 
diverges, then eq. ( |5.5| ) implies that d — 2 + r] < 0, so that the critical connected 2-spin 
correlation function, eq. (|5.2|), diverges with the distance r between the spins, rather than 



■^As usual, logarithmic factors could, in principle, be present, in particular, at the upper critical dimen- 
sionality; they would not change the value of the algebraic exponents. 
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having its usual decay to zero as r ^ oo. Of course, this behaviour could never happen for 
physical temperatures, but it is forced on us at a complex-temperature critical point where 
(i) M diverges, provided that (ii) the correlation length also diverges at this point and, (ii) 
the usual scaling assumption for the 2-spin correlation function is valid. 

6 Behaviour of the Specific Heat 
6.1 General 



By re-expressing the exact solution |00 for the specific heat C for the Ising model on the 



,-^^K-^C = -,A^ + ^'""' + !,l''":^'*- /a^-<) - 5(l^#iE(*,) (6.1.1) 



triangular lattice in the FM phase in terms of the variable u, we obtain 

where K{k) = J^^^ de[l - k'^ sin^ 0]-V2 and E{k) = J^^^ de[l - ^ sin^ ^] V2 are the complete 
elliptic integrals of the first and second kinds, respectively, depending on the the elliptic 
modulus k, and A;< was given in eq. ( |3.4.1| ). 



6.2 Vicinity of u = —1/3 

To consider the approach to the point u = —1/3 from within the complex extension of the 
FM phase, we use eq. (|3.4.2| ) with e ^ and hence ( p.4.3|) . We let /c< = in/n', where n' is 



the complementary elliptic modulus satisfying k^ + k'^ = 1. It follows that as n ^ —1/3) 
K = e/(4 + e) ^ 0. Using the identity [|| 



K^in/n') = k'K{k) (6.2.1) 

we find that the term involving K{k^) in eq. (|6.1.1|) approaches a constant times (l+3'u)~^'^. 
In the second term, one factor of (1 + 3'u)^^/^ comes from the /c< while another comes from 
the elliptic integral E{k^), so that this term diverges like (1 + 3n)~^. This is the leading 
divergence in C, so we thus obtain the exact result that as m — > —1/3 from within the 
complex FM phase, the critical exponent for the specific heat is 

a', = 1 (6.2.2) 



''Houtappel's expressions for the internal energy and specific heat, eqs. (108) and (109), respectively, 
in Rcf. |l^, are incorrect if one uses the integrals ei(/3) and 62 (/3) as he defines them, with the range of 
integration from (j) = to 2-k. If, instead, one takes the range of integration from cj) — Q to (p = 7r/2, so that 
the integrals are just the usual elliptic integrals K{^/]3) and E{^/]3), then his eqs. (108) and (109) become 
correct. 
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To our knowledge, this is the first time an algebraic power has been found for the specific 
heat critical exponent at a singular point in a 2D Ising model. For the critical amplitude, 
we calculate (taking = in eq. (|3.4.2|) ) 

k^'K-'C ^ -M^|l + 3u\-' (6.2.3) 

n 

The infinite set of K values corresponding to the point m = — 1/3 is 

/f = iln3-P" + '>"^ (6.2.4) 

4 4 ' ^ 

where n E Z. 

6.3 Vicinity of u = —1 

We first observe that as m ^ — 1 from within the complex FM phase, /c< ^ — 1. It follows 



that in this case the term involving E(k^) in eq. ( 6.1.1 ) is finite while the term involving 



K{k^) diverges logarithmically, so that at this point Ug = —1, 

a^ = {log div) (6.3.1) 

Thus, the divergence in the specific heat at u = — 1 on the triangular lattice is of the same 
logarithmic type that it is |^, |T^ at -u = — 1 on the square lattice, and the same as it is at the 
respective physical critical points on both the square and triangular (as well as honeycomb) 
lattices. For the critical amplitude, using the Taylor series expansion of fc<, as a function of 

-1, 

fc< = -1 - 2-='(l + uf + 0((1 + uf) (6.3.2) 



we calculate 

12i 
ks^K~^C^ — ln|l + M| as u ^ -1 (6.3.3) 

n 

The infinite set of K values corresponding to this point was given in eq. (3.1.8) of paper I: 
K = —in/A + niix /2 where n E Z . 

One can also consider the approach to m = —1 from within the complex-temperature 
extension of the symmetric, PM phase (from the upper left or lower left in Fig. 1(a)). Using 
the exact expression for the specific heat applicable in the symmetric phase [jl3|, Q , we find 
the same logarithmic divergence in C, so that the corresponding exponent is 

as = {log div) (6.3.4) 
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7 Behaviour of x ^it the Boundary of the PM Phase 



For the square lattice, we have proved g, |T^ that the susceptibihty has at most finite 
singularities on the boundary of the complex FM phase, aside from the physical critical 
point and the singular point at m = — 1. We expect a similar result to hold for the triangular 
lattice, although to show this with complete rigour, it would be desirable to perform an 
analysis of the asymptotic behaviour of the general connected 2-spin correlation function 
< cro,o<^m,n >, as r = [m? + n^)^/^ — > oo, which has not, to our knowledge, yet been done. 

We have carried out a study of the high-temperature expansions for x on the triangular 
lattice in order to search for complex-temperature singularities. Since these singularities 
are presumably finite, we used the method of differential approximants, which is capable 
of representing such a finite singularity in a function even in the presence of an analytic 
background term. The high-temperature series for x on the triangular lattice was calculated 



in 1971 to 0{v ) ||3^. We have checked and found that apparently it has not been calculated 



to higher order since p6| . Mainly because of the shortness of the series, our study did not 



yield a definite value for the exponent 7^ at the points f = ±i corresponding to the point 
M = — 1 as approached from within the complex PM phase. This might be possible with a 
substantially longer high-temperature series. 

8 Conclusions 

In this paper we have investigated the complex-temperature singularities of the (zero-field) 
susceptibility x for the Ising model on the triangular lattice. As part of this study, from 
an analysis of the low-temperature series using differential approximants, we have obtained 
evidence which agrees with and strengthens the earlier conclusion by Guttmann [|l[ from 
ratio and dlog Fade methods that x diverges like |1 — u/ue\~'^'' at the endpoint Ue = —1/3 
with exponent 7g = 5/4 and have obtained a somewhat more precise determination of the 
critical amplitude, again in agreement with Ref. |l[]. We also discussed how this behaviour 
of X is related to the unusual feature that the spontaneous magnetisation diverges at this 
point, and some implications of this. From exact results, we extracted the specific heat 
critical exponents at m = —1/3 and u = —1, which are, respectively, a'^ = 1 and a^ = (log 
divergence). Finally, we have commented on non-analyticities of x which could occur as one 
approaches the boundary of the symmetric phase from within that phase. 

One of us (RS) would like to thank Prof. David Gaunt and Prof. Tony Guttmann for 
information about the current status of the series expansions for the triangular lattice. RS 
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would also like to thank Prof. Tony Guttmann for kindly informing us of several of his early 
works on complex-temperature singularities and for discussions of these works. This research 
was supported in part by the NSF grant PHY-93-09888. 
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Figure Caption 

Fig. 1. Phases and associated boundaries in the complex variables (a) u, (b) z, and (c) 
V, as defined in eqs. (|2.2| )-( ^73|) . In Fig. 1(a), the locus of points is given by the union of eqs. 
( p.ll| ) and ( p.l2| ) in the text. The intersections of the circle with the real u axis occur at 
Uc = 1/3 and Ug = —1, the latter of which is also its intersection with the semi-infinite line 
segment (|2.12|) extending from the point u = —1/3, denoted Ue in the text, to u = — oo along 



the negative real axis. In Fig. 1(b), the intersection of the outer boundary of the complex 
FM phase with the real z axis consists of the two points ±Zc = ±l/-\/3- The two vertical line 
segments extend, respectively, from ±i/^/3 to ±ioo along the positive (negative) imaginary 
axis. In Fig. 1(c), the intersections of the curves with the real v axis occur, from left to 
right, at f = —1, v = Vc = 2 — \/3, and v = v~^ = 2 + -\/3. The endpoints of the arcs occur 
at V = exp(±27r/3). 
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